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ABSTRACT 



Aims. We determine the cumulative spectral energy distribution (SED) of local active galactic nuclei (AGN) in the 
3-300 keV band and compare it with the spectrum of the cosmic X-ray background (CXB) in order to test the widely 
accepted paradigm that the CXB is a superposition of AGN and to place constraints on AGN evolution. 
Methods. We performed a stacking analysis of the hard X-ray spectra of AGN detected in two recent all-sky surveys, 
performed by the IBIS /ISGRI instrument aboard INTEGRAL and by the PGA instrument aboard RXTE, taking into 
account the space densities of AGN with different luminosities and absorption column densities. 

Results. We derived the collective SED of local AGN in the 3-300 keV energy band. Those AGN with luminosities below 
1043-5 gj.g g-1 (■i7_5o keV) provide the main contribution to the local volume hard X-ray emissivity, at least 5 times 
more than more luminous objects. The cumulative spectrum exhibits (although with marginal significance) a cutofi^ at 
energies above ~ 100-200 keV and is consistent with the CXB spectrum if AGN evolve over cosmic time in such a way 
that the SED of their collective high-energy emission has a constant shape and the relative fraction of obscured AGN 
remains nearly constant, while the AGN luminosity density undergoes strong evolution between z ~ 1 and z = 0, a. 
scenario broadly consistent with results from recent deep X-ray surveys. 

Conclusions. The first direct comparison between the collective hard X-ray SED of local AGN and the CXB spectrum 
demonstrates that the popular concept of the CXB being a superposition of AGN is generally correct. By repeating 
this test using improved AGN statistics from current and future hard X-ray surveys, it should be possible to tighten 
the constraints on the cosmic history of black hole growth. 
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1. Introduction 

It is widely believed that the bulk of the cosmic X-ray back- 
ground (CXB) consists of emission from all active galac- 
tic n uclei (AGN) in the visible Universe (jSetti fc Woltieil 
. Il989t ). Therefore, the CXB represents a unique integral 
' record of the history of growth of massive black holes since 
' the early epochs until the present time. The popularity of 
this concept rests on the absence of significant distortions in 
' the thermal spectrum of the cosmic microwave background, 
, which rules out that a hoi : intergalactic medi um is a major 
• contributor to the CXB ijWright et al.lll994[) . and also on 
the fact that most of the background emission at energies 
below sev eral keV has been directly resolved into individual 
AGN (see lBrandt fc Hasingeij|2005l for a recent review). 

However, the peak of the CXB spectral energy distri- 
buti on (SED) is situated near 30 keV (e.g. iGruber et al.l 
|199S) . whe re ^ 99% of the bac kground emission remains 
unresolved ([Krivonos et al.l[2007[ ) . Fo r this reason, cur rently 
popular CXB synthesis models (e.g. iGilh et aI]|2007D com- 
pare versatile AGN statistics (source counts, luminosity 
functions at different redshifts, distribution of absorption 
column densities, etc.) collected by X-ray surveys at ener- 
gies below 8 keV with the CXB spectrum, under certain 
assumptions about AGN SEDs in the hard X-ray and soft 



gamma-ray bands (from 10 to several hundred keV). This 
allows one to test the overall paradigm of the CXB being 
a superposition of AGN and possibly to get an idea about 
still missing (e.g. heavily obscured) populations of AGN 
and other types of extragalactic objects. 

Recent observations have provided a number of con- 
straints for this picture. In particular, it has been estab- 
lished how the AGN luminosity function defined at rest- 
frame energies below 8 keV evolves with redshift, at least 
out to z ~ 2, and there are now constraints on the frac- 
tion of ob scured AGN as a function of luminosity a nd red- 
shift (e.g. lUeda et al.l 120031: iBarger et all 120051 : iTozzi et all 
i200§). Most of these data were obtained in deep, pencil- 
beam extragalactic surveys conducted by the Chandra and 
XMM-Newton observatories. 

Valuable complementary statistics have recently 
been provided by shallow hard X-ray surveys of the 
wh ole sky performed b y RX TE (iRevniytsev et all 120041: 
Sa zonov fc RevnivtsevI 2004f). Swift (iMarkwardt et al 
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200JI) , INTEGRAL (iBassani e t al.' 2006'; 'BeckmanneLaL 
2006t iKrivonos et all 120071: .Sazonov et al.. ,20071 ) and 
HEAD 1 (|Shinozakill2006ir Most importantly, these surveys 
have made it possible for the first time to reliably measure 
the distribution of AGN absorption column densities well 
into the Conipton thick regime (A'h ^ 10^^ cm~^), albeit 
only at low redshifts (z < 0.1). In particular, it was found 
that the observed ratio of obscured (A'h > 10^^ cm^^) to 



unobscured (iVn < 10 cm""^) AGN drops from about 2:1 
at hard X-ray (17-60 keV) lumino sities ;$ lO^'^'^ erg to 
abo ut 1:3 at higher luminosities (jSazonov ct al. 200 7t see 
also ISazonov fc RevnivtsevI 12004 iMarkwardt et al.i 120051 : 
IShinozakill20'06l^ 

As was mentioned above, previous CXB studies used 
some fiducial intrinsic AGN SED going up to several hun- 
dred keV, which was usually assumed to be similar to the 
few measured spectra of brightest nearby Seyfe rt galaxies 
(e.g. iJourdain et al.|[l99^ IZdziarski et aLlfTQQSf) . However, 
it has never been demonstrated that the cumulative hard 
X-ray SED of all AGN residing in a given volume of 
the Universe is compatible with the CXB spectrum. It is 
only now that such crucial comparison can be made for 
the first time, albeit only for the low-redshift AGN pop- 
ulation, using the all-sky hard X-ra y surveys mentioned 
above. The IBIS/ISG RI instrument lubertini et al.l l2003D 
aboard INTEGRAL (jWinkler et all l2003f ) is particularlv 
suitable for this purpose since it is effectively sensitive up to 
300 keV, i.e. well into the energy range where the collective 
AGN SED is expected to have a cutoff if our understanding 
of the CXB origin is correct. 

The purpose of the present work is to estimate the cu- 
mulative hard X-ray SED of local {z <, 0.1) AGN. For this 
purpose we stack the spectra of AGN detected in two recent 
all-sky surveys, performed by the IBIS/ISGRI instrument 
aboard INTEGRAL and by the PCA instrument aboard 
RXTE, taking into account the space densities of AGN with 
different luminosities By comparing the derived cumu- 
lative SED of local AGN, which spans two decades in energy 
(3-300 keV), with the CXB spectrum we obtain constraints 
on the evolution of AGN over cosmic time (fj3]). A cosmol- 
ogy with r^in — 0.3, r^A — 0.7, and Hq = 75 km s^^ Mpc^^ 
is adopted throughout the paper. All quoted uncertainties 
are la unless noted otherwise. 



2. Analysis 

We recently used the (mostly) serendipitou s all-sky survey 
conducted by the IBIS/ISGRI instrument (|Krivonos et al.l 
dOOl) to obtain a samp le of nearby AGN d etected in the 
17-60 keV energy band (jSazonov et al.ll2007D . This sample 
is well suited for estimating the cumulative hard X-ray SED 
of local AGN. In what follows we first describe our stacking 
spectral analysis performed on this AGN set f H2.ip . We then 
report on our similar analysis carried out at lower energies 
(3-20 keV) using the RXTE Slew Survey Finally we 

put both sets of results together to obtain a composite AGN 
SED covering the energy range 3-300 keV {i 



2.1. INTEGRAL sample 



The all-sky survey reported by iKrivonos et al.l (|2007T ) is 
based on INTEGRAL/IBIS/ISGRI observations performed 
during 2002-2006, including a special series of observations 
of "empty" extragalactic fields. It greatly improves on pre- 
vious hard X-ray surveys in terms of angular resolution 
(~ 12') and sensitivity. For 80% of the sky a source detec- 
tion limit of 5 mCrab « 7x 10"^^ erg s~^ cm"^ (17-60 keV) 
or better is achieved. A total of 403 sources were found, with 

at most 1-2 of them being spurious. 

Ne arly two thirds of the sources in the IKrivonos et al.l 
(|2007l ) catalog reside in the Galaxy, while the rest are 



confirmed or suspected AGN. We previously used the 
INTEGRAL catalog to construct the hard X-ray (17- 
60 keV) luminosity function of nearby AGN and to 
study their distribu tion in absorption column density 
(jSazonov et al.|[2007t ) . As in that work, we now use the sub- 
sample of AGN located outside the Galactic plane region 
{\b\ > 5°) to avoid problems with identification incomplete- 
ness at low Galactic latitudes; our AGN sample is highly 
complete (~ 95%) at |5| > 5°. 

In fact the sample of AGN used in th e present work 
is shghtly different from that presented by ISazonov et al.l 
(20 (33) but p r ecisely corresponds to the catalog of 
iKrivonos et al.l (j2007f ). which included a few additional 
INTEGRAL observations and follow-up identifications that 
had become available between those two publications. 
Specifically, our current sample is composed of 76, rather 
than 74 AGN detected on the IBIS/ISGRI average sky map 
at |6| > 5°. Information on the two additional sources is 
presented in Table [11 which should be reg arded as a contin- 
uation of Table 1 iii lSazonov et al.l (j2007j ). Among these 76 
AGN, there are 8 blazars and 68 nearby {z ^ 0.1) Seyfert 
galaxies (although a few of them still lack an exact optical 
classification). Below we focus on these non-blazar AGN, 
although a short notice is made in ii2.4l with respect to the 
blazar contribution to the cumulative AGN spectrum. 

2.1.1. Reconstruction of multiband sky images 

To obtain the hard X-ray spectra of the INTEGRAL AGN 
we essentially repeated our analysis of IBIS/ISGRI data, 
previously performe d in the 17-60 keV energy band (see 
IKrivonos et cd]l2007l for details), in 7 narrow channels: 17- 
26, 26-38, 38-57, 57-86, 86-129, 129-194 and 194-290 keV. 
This resulted in time-averaged maps of the whole sky in 
these energy bands, which were then used to measure source 
spectral fluxes. In this analysis, the source posit i ons w ere 
fixed at the values determined bv IKrivonos et al.l (j2007l ) in 
the 17-60 keV band. 

There are two imp ortant differences with respect to the 
IKrivonos et al.l (j2007[ ) study. First, for the present work 
we updated our data set by adding all IBIS/ISGRI data 
that became available to us since the last publication. This 
amounted to 26 Ms of cleaned and deadtime-corrected data 
in addition to 33 Ms available before. Secondly, in recon- 
structing clean sky images from raw IBIS/ISGRI shadow- 
grams in a given energy channel we used only those sources 
detectable in that band rather than the entire source cata- 
log. This essentially eliminates iterative removal of sources 
from images, unnecessary in the high-energy channels (129- 
194 and 194-290 keV), where only 27 and 8 sources, respec- 
tively, are detected (> 5(t) on the whole sky, and reduces 
the error in measured fluxes. 

It is important to emphasize that despite the addition 
of the most recent IBIS/ISGRI data we co ntinued to base 
our a nalysis on the original AGN catalog (jSazonov et al.l 
|2007| plus the two so urces in Table [B and t he corresponding 
exposure map from IKrivonos et alj (j2007D . Thus, the new 
data, added in an attempt to improve the quality of the 
source spectra, do not affect the statistical properties of 
our AGN sample. 

When carrying out a stacking analysis like the one re- 
ported here one should always worry that some systematic 
uncertainties associated with the reconstruction of individ- 
ual source fluxes may be greatly enhanced by adding up 



Table 1. AGN detected on the average IBIS/ISGRI map at \b\ > 5° in addition to the lSazonov et aP (|2007l ) sample 



Object 


Class'^ 


Rof." 


z 


D" 






Nn 


Rcf." 










Mpc 


10"^^ erg s"^ cm"^ 


erg ^ 


10^^ cm"^ 




SWIFT J0601.9-8636=ESO 005-G004 


S2 


1 


0.0062 


22.4 


2.51 ± 0.46 


42.18 


~ 100 


1 


ICR J14561-3738=ESO 386-G034 


S2 


2 


0.0246 




1.40 ± 0.26 


43.23 


> 100 


2 



^ Optical AGN class: S2 - Scyfcrt 2 galaxy. 

^ Reference for the optical classification: (1) IMorelli et all f200^. (2 ) [Masetti et al.| (|2007f) . 
Distance according to the Nearby Galaxies Catalogue flTullvHigSST) . 

Observed (uncorrected for absorption) luminosities in the 1 7—60 keV energy b and. 

® Quoted absorption column density A'^h is adopted from: fl) lUeda et al.| | |2007D , f 2) ISazonov 



the fluxes of many sources. Therefore, to verify that our 
results do not suffer from systematic uncertainties related 
to image reconstruction we made a number of simulations. 

In particular, we have extracted spectra from 3,000 
"empty" positions on the sky chosen in accordance with the 
INTEGRAL exposure map, i.e. the number density of these 
"zero-flux" sources was made higher in regions of large ex- 
posure than in regions of small exposure. When we carried 
out simulations including in the iterative source removal 
procedure all the sources from the INTEGRAL catalog, 
irrespective of the significance of their detection in the con- 
sidered energy channel, the mean of the distribution of the 
fluxes of the simulated sources proved to be ^ O.OS-O.lo- 
above zero (where a is the typical uncertainty of source 
flux measurement), which is unacceptable for our stacking 
analysis. This result is actually antici pated for the ernploye d 
algorithm of image reconstruction (jKrivonos et al.ir2007f) . 
However, when we repeated the simulations in exactly the 
same manner as we carried out our real stacking analysis, 
i.e. including into our iterative source removal procedure 
only those INTEGRAL sources that are signiflcantly de- 
tected in the considered energy channel, the distribution of 
the fluxes of the simulated sources proved to be consistent 
with the normal one with zero mean and the dispersion 
expected from Poisson statistics. 

This proves that a stacking analysis of as many as 3,000 
sources does not suffer from any systematic uncertainties 
in our image reconstruction procedure. 

2.1.2. Stacked spectra 

The individual AGN spectra derived from the reconstructed 
IBIS/ISGRI multiband sky images will be discussed in de- 
tail elsewhere. Here we are only interested in the average 
properties of these spectra. 

We first separately stacked AGN spectra within three 
groups defined by the source signal-to-noise ratio in the 
17-60 keV energy band: > 30 (5 sources), between 15 and 
30 (6 sources) and < 15 (57 sources). This is done mainly 
to prevent the stacked spectrum from being dominated by 
one or two bright sources. The resulting spectra are shown 
in Fig. [TJ One can see that all three spectra are similar to 
each other. They exhibit some rollover between ~ 40 and 
~ 200 keV and can be approximately described by a power 
law with a high-energy cutoff: dN/dE oc exp {—E/Ef), 
with r ~ 1.65 and Ef 250 keV (see the solid lines in 
Fig- El)- This model provides an adequate approximation of 
the observed spectra, but given the limited energy coverage 
(17-290 keV), other functional forms (e.g. a broken power 
law) would flt the INTEGRAL data equally well. 

The fact that the averaged spectra of local AGN pre- 
sented in Fig.[l]have a maximum at energies ~ 80-100 keV, 
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Fig. 1. Average IBIS/ISGRI spectra of non-blazar AGN 
divided into three groups according to their detection sig- 
niflcance in the 17-60 keV energy band. The spectral nor- 
malizations are arbitrary. Upper limits are la. The solid 
lines show crude approximations by the cutoff power law 
model with F = 1.65 and Ef = 250 keV. 

which can be explained by the presence of a high-energy 
cutoff with Ef ^ 250 keV, implies that if the same hard X- 
ray emission were coming from redshifts z ~ 1-2 its spec- 
trum observed a.t z = would have a peak near ~ 30 keV, 
similar to the cosmic X-ray background. However, these 
stacked spectra cannot yet be regarded as characteristic of 
the collective hard X-ray emission of the local AGN popu- 
lation, since neither AGN space densities nor luminosities 
were taken into account when constructing them. 

To really obtain the cumulative SED of local AGN, the 
following weighted stacking should be done: 



S^ 



Vniax,7 
3 J 



j2 ^^D^p^^j 



(1) 



Here Si is the cumulative volume emissivity (measured in 
units of erg s~^ Mpc~^) in energy channel i (from 1 to 
7), Li J is the observed luminosity of j'th AGN in chan- 
nel z, Fi^j is the measured flux from j'th AGN in channel 



i, Dj is the source luminosity distance and Knaxj" is the 
volume of space within which an AGN with observed lu- 
minosity ijixj- in the 17-60 keV band could be detected in 
the INTEGRAL survey. We adopted the source distances 
{Dj) and signal-to-noise ratios (used to derive the maxi- 
mum volu mes T^nax.? given the INTEGRAL sky exposure 
map) from lSazonov et al] ()2007f ). where the hard X-ray lu- 
minosity function of local AGN was derived in a similar 
way, and also from Table [1] The conversion of detector 
counts to photon fluxes in the 7 spectral bands was cali- 
brated using IBIS/ISGRI measurements of the Crab and 
assuming that the Crab spectrum is given by dN/dE — 
10(i?/l keV)^^-^ phot keV~^ s~^ cm~^. In calculating the 
Vma^j volumes we also assumed the observed spectra to be 
Crab-like. Although in general this may be a poor assump- 
tion for heavily obscured (logTVn » 10^^ cm~^) AGN, it 
is unlikely to introduce a significant error in our case, be- 
cause we use a relatively narrow energy band (17-60 keV) 
for source detection and since such strongly absorbed AGN 
provide a small contribution to the resulting cumulative 
SED. 

There are two types of uncertainty associated with the 
summed spectrum 5"^. One results from the errors SFij in 
the measured source spectral fluxes (and the associated un- 
certainties, SLij = AnD^SFi j , in the spectral luminosities) 
and is given by: 



E 



max, J 



(2) 



Another uncertainty is associated with the finite size of our 
AGN sample: 



5S, 



i,2 



\ 



L, 



max, J 



(3) 



Therefore, these uncertainties are associated with the 
Poisson distribution of detector counts and the number of 
sources, respectively. Since the first uncertainty is directly 
linked to the survey's sensitivity in different energy bands, 
its relative amplitude is much larger for the high-energy 
channels than for the low-energy ones. In contrast, the sec- 
ond uncertainty is nearly independent of the energy band 
{SSi^2/Si ~ const) and affects the normalization of the cu- 
mulative SED rather than its shape. 

In Figs. [H and [3] we show two cumulative AGN SEDs 
obtained using equations ([I])-® - In the first case the stack- 
ing was done for the AGN (40 in total) with luminosities 
between L^in = 10''^ and L^ax = lO*^-^ erg s'^ (17- 
60 keV). In the other case, Ly^in = lO'''^'^ erg s~^ and 
imax = oo (27 AGN). This division approximately corre- 
sponds to the low-luminosity and high-luminosity branches 
of the luminosity function of nearby AGN, which are char- 
acterized bY_differeiit ratios of obscured and unobscured 
sources (jSazonov et al.l [20071) and so may represent some- 
what different populations of AGN. Note that one source in 
our sample - the Seyfert galaxy NGC 4395 - has luminosity 
less than 10'*^ erg s~^ (Lhx ~ 40.4) and was not used in the 
stacking analysis. 
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Fig. 2. Cumulative 3-300 keV SED of local AGN with lumi- 
nosities (observed in the 17-60 keV band or absorption cor- 
rected in the 3-20 keV band) 10^^ < Lhx < lO''^'^ erg s~\ 
obtained with INTEGRAL/IBIS/ISGRI and RXTE/PCA. 
The error bars are la uncertainties of the first type [equa- 
tion while the shaded regions show the combined la 
uncertainties of the first and second type [equation ([3])]. 
The solid line is the best-fit model by a sum of absorbed 
and unabsorbed power laws with a high-energy exponential 
cutoff [equation (|4]), Table [4]. The dashed lines indicate the 
contributions of unobscured AGN, and obscured AGN with 
A^H < 10^'' and > lO^'' cm'^. The RXTE points have been 
mutliplied by a factor of 1.1 and the INTEGRAL points 
divided by the same factor to correct for the effect of the 
local large-scale structure (see t j2.2.ip . 



2.2. RXTE/PCA sample 

We have thus obtained the cumulative SEDs of local low- 
luminosity and high-luminosity AGN in the energy range 
17-290 keV. There is an obvious need to continue these 
spectra to lower energies where intrinsic X-ray absorption 
in AGN should come into play, producing a low-energy cut- 
off in the spectra of obscured objects. There are two possi- 
ble ways to obtain this low-energy part of the cumulative 
SED. One is to search the literature and data archives for 
X-ray spectra of the INTEGRAL AGN. We decided not 
to follow this approach for two reasons. First, published 
or archival X-ray observations of sufficient quality do not 
exist for all of our objects. Second, because AGN are inher- 
ently variable, adding lower-energy spectral fiuxes to the 
IBIS/ISGRI hard X-ray spectra could introduce a bias due 
to the fact that these particular AGN tended to be detected 
by IBIS/ISGRI as they were bright relative to their typical 
flux levels (averaged over many years), while a different X- 
ray telescope will find the same source equ ally likely either 
below or above its average flux level fsee IShinozakil |2006| 
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Fig. 3. Cumulative 3-300 keV SED of local AGN with 
Lhx > lO**^'^ erg s^^; see Fig. [5] for explanations. The up- 
per limit in the highest energy channel is la. We note the 
apparent mismatch between the INTEGRAL and RXTE 
parts of the SED, which is statistically insignificant and 
may be attributable in part to effects of the local large- 
scale structure, which we cannot accurately estimate and 
thus have not corrected for in this high-luminosity case (see 
in giXl]and ^O]) . 

for a detailed discussion of this effect). As a result, the low- 
energy part of the cumulative SED derived by this method 
would somewhat underestimate the true SED of the local 
AGN population. 

Another feasible approach consists of using an indepen- 
dent survey of nearby AGN performed at the energies of 
our current interest, i.e. just below 20 keV. Such a survey 
of the whole sky in the energy band 3-20 keV was recently 
performed with the PCA instrument aboard the RXTE ob- 
servatory during its slews betwee n pointed observa.tions - 
the RXTE Slew Survey, or XSS (|Revnivtsev et al.ll2004n . 
In this survey a flux limit of 1.8 x 10~^^ erg s~^ cm~^ 
(3-20 keV) or better was achieved for 80% of the sky at 
\b\ > 10°. Therefore, the XSS is substantially more sen- 
sitive than the IBIS/ISGRI survey with respect to unob- 
scured AGN (e.g. Seyfert 1 galaxies) but is less sensitive 
with respect to strongly obscured (TVh 10^'^ cm~^) ones. 

The high Galactic latitude (|5| > 10°) part of the XSS 
was previously used to obtain a catalog of AGN, to con- 
struct the 3-20 keV luminosity function of nearby AGN 
and t o study their distribution of a bsorption column den- 
sities (jSazonov fc RevnivtsevI [200^ . quite similarly to the 
analysis of the INTEGRAL survey we carried out later. 
Here we use the XSS to obtain low-energy intensities for 
our cumulative hard X-ray SEDs of local AGN. 

Since the original XSS publications several changes have 
taken place that need to be taken into account here. First, 



originally a substantial fraction (~ 30%) of XSS sources 
were unidentified, mainly because of their poor localization 
(uncertainty up to 1°). Over the passed three years many 
of these sources h ave been identif i ed, thanks to dedicated 
follow-up efforts (iBikmaev et al.l l2006t iRevnivtsev et al.l 
120061 : iLandi et al.ll2007l ) or due to detection and improved 
localization of some of these sources by INTEGRAL and/or 
Swift. As a result, the XSS sample of AGN detected (> Aa) 
at \b\ > 10° is now composed of 103, rather than 95 sources, 
including 84 non-blazar AGN. The new non-blazar AGN 
are listed in Table [2l which should be considered a con- 
tinuation of Table 1 in ISazonov fc RevnivtsevI ()2004[ ). At 
the same time the number of unidentified XSS sources (de- 
tected at > 4(t) has decreased to 16 objects. Therefore, the 
XSS AGN sample is now at least 82% complet43, so incom- 
pleteness is significantly less of a problem than before. 

Secondly, through careful testing of our methods devel- 
oped for analyzing data recorded in the slew mode of RXTE 
observations we found that the conver sion factors from 
PCA counts to source photons quoted in IRevnivtsev et all 
()2004D had been underestimated by a factor of 1/0.7, be- 
cause we had not taken into account that the PCA field of 
view is rapidly moving across sources during RXTE slews 
and used a slow movement approximation of the PCA re- 
sponse. This implies that, to a first approximation, the 
AGN luminosities quoted in ISazonov fc RevnivtsevI (|2004f ) 
should be revised upwards by a factor of 1/0.7 (the same 
correction factor applies to Table [2] presented here) . 

There are a few other minor modifications that have 
been implemented, including the correction of a misprint in 
ISazonov fc RevnivtsevI (I2004D for the redshift of NGC 7582 
from 0.053 to 0.0053 and accordingly revising its luminosity. 

We plan to release an updated version of the XSS cata- 
log in the near future, including the modifications described 
above and a number of additional source identifications. 

It is important to note that although the RXTE and 
INTEGRAL AGN samples are comparable in size, they 
overlap by just ~ 30%, i.e. these samples are mostly in- 
dependent of each other. Although this result may seem 
unexpected, it can readily be explained by a combination 
of circumstances: 1) the XSS catalog is defined at |6| > 10°, 
whereas the INTEGRAL AGN sample used here is defined 
at |6| > 5° and the strip of the sky between |6| = 5° and 
|6| = 10° has been well covered by INTEGRAL observa- 
tions, 2) the distribution of exposure over the sky is sig- 
nificantly different for the two surveys, 3) due to its lower 
energy band the XSS is biased toward detecting unobscured 
and weakly obscured AGN and against detecting Compton 
thick AGN, 4) AGN are inherently variable and the two 
surveys were conducted at significantly different epochs. 

We used the updated XSS AGN sample to obtain the 
3-8 and 8-20 keV data points for our cumulative hard X- 
ray SEDs of local low-luminosity and high-luminosity AGN 
(Figs, m and 13]) following the same 1/Vmax weighted sum- 
ming method as we used above f H2.ip to obtain the high- 
energy parts of the SEDs based on IBIS/ISGRI data. We 
note that the XSS sample includes 6 non-blazar AGN lo- 
cated at z ~ 0.15-0.3. Nonetheless, since the vast majority 
of the sample are truly local sources (z ^ 0.1) and the sev- 
eral more distant AGN contribute just a few per cent to 
the resulting cumulative SEDs, we included these objects 



^ Here it is also taken into account that ~ 6 additional AGN 
are possibly confused with other XSS sources. 



Table 2. Non-blazar AGN detected at |5| > 10° during the RXTE Slew Survey in addition to the lSazonov fc RevnivtsevI 
(|2004f) sample 



XSS object 


Name 


Class^ 


Ref." 


3-8 keV 


8-20 keV 


z 




-20*^ 




Ref* 


(J2000.0) 








ent s^^ 


ent s^^ 




erg s 




10^^ em-2 




05054-2348 


2MASX J05054575 — 2351139 


S2 


1 


0.69 ± 0.19 


0.92 ± 0.23 


0.0350 


43.77 


43.83 


6 


1 


12303-4232 


IRAS F12295-4201 


SI. 5 


2 


0.48 ± 0.09 


0.29 ± 0.11 


0.1000 


44.37 


44.37 


< 1 


1 


12389-1614 


2MASX J12390630-1610472 


S2 


3 


0.93 ± 0.11 


0.51 ± 0.13 


0.0367 


43.75 


43.77 


2 


2 


15076-4257 


2MASX J15080462-4244452 


SI 


4 


0.73 ± 0.05 


0.24 ± 0.06 


0.0565 


43.96 


43.96 


< 1 


3 


18236-5616 


IC 4709 


S2 


5 


0.32 ± 0.09 


0.51 ± 0.12 


0.0169 


42.86 


42.98 


12 


1 


19459+4508 


2MASX J19471938+4449425 


S2 


1 


0.38 ± 0.10 


0.38 ± 0.12 


0.0532 


43.84 


43.95 


11 


2 


21354-2720 


IRAS F21318-2739 


SI. 5 


1 


0.41 ± 0.11 


0.32 ± 0.13 


0.0670 


43.99 


43.99 


< 1 


4 



^ Optical AGN class: SI. 5 - Seyfert 1.5 galaxy, S2 — Scyfcrt 2 galaxy. 

^ Reference for the optical classifi cation: ( 1) Bikmacv ct al] j2006D . (2) ILandi et all ) |2007D , (3) IMasetti eFall (|2Q06aD , (4) a Seyfert 1 nucleus is 
suggested by the 6dF spectrum, (5) IMasetti et al.. (ii2006b,) . 

^ Left and right columns give the observed and absorption corrected lumino sities in the 3—20 keV ba nd, respectively. 

Quoted absorption column density value is adopted from or based on: fl) iRevnivtsev et al.l (]2006r ).f2) ISazonov et al.l ll2005h . (3) ROSAT data, (4) 
Swift data. 



in our stacking analysis as if they were local, i.e. ignoring 
any evolution of AGN between z = 0.3 and z = 0. 

There are three differences with respect to our analysis 
of the IBIS/ISGRI data, all caused by the strong effect of 
intrinsic absorption on AGN spectra below 20 keV. First, 
in determining the Knax values for the XSS sources, we 
took into ac count their absorption column densities given 
in Table 1 in lSazonov fc RevnivtsevI (|2004 and in Table [2] 
presented here. That is Knax is now defined as the volume 
of space in which the XSS would detect an AGN with its 
observed 3-20 keV luminosity and column density iVn (not 
just the luminosity as in the INTEGRAL case). Second, 
the TVh values were taken into account also to calculate 
the observed and intrinsic (corrected for the absorption) 
luminosities of the XSS AGN in the subbands 3-8 and 8- 
20 keV from the detec tor counts measured in these c hannels 
(as given in Table 1 in lSazonov fc RevnivtsevI (|2004l l and in 
Table [5] presented here). In doing this conversion the source 
spectrum was assumed to have an intrinsic slope of 1.8, 
although the result only weakly depends on this assumption 
(jRevnivtsev et aLll2004l ) . And finally, we used the same two 
luminosity ranges: 1) Lmin = 10"'-^, imax = lO^"^-^ erg s~^ 
(38 AGN) and 2) L^in = 10^^ \ ^max = oo (46 AGN), 
as for our stacking analysis of INTEGRAL spectra, but for 
RXTE these ranges are defined for intrinsic (i.e. absorption- 
corrected) luminosities in the 3-20 keV energy band. This 
ensures that the stacking of both INTEGRAL and RXTE 
spectra is done for almost identical populations of AGN, 
since the unabsorbed luminosities of AGN in the 3-20 keV 
and 17-60 keV ban ds are very similar (see e.g. Fig. 7 in 
ISazonov et "alll2007[ ). 

2.2.1. Influence of the local large-scale structure 

Both the INTEGRAL and RXTE all-sky surveys are char- 
acterized by substantially nonuniform coverage of the sky. 
Furthermore, the sky exposure maps of these surveys are 
quite different. Because both surveys effectively picked out 
AGN within ^ 100-200 Mpc from us and the distribution of 
matter in the local Universe is inhomogeneous, effects of the 
local large-scale structure (LSS) may come into play when 
comparing AGN space densities or any quantities depend- 
ing thereof (luminosity functions, cumulative SEDs etc.) 
estimated using one survey and the other. It is straightfor- 
ward to estimate these effects if we assume that AGN are 
distributed in the local Universe approximately as normal 



galaxies. iKrivonos et al.l (|2007l ) have actually demonstrated 
that the AGN detected by INTEGRAL do follow quite well 
the spatial distribution of galaxies. 

To assess the effect of the local LSS on the cumulative 
SEDs of low- and high-luminosity AGN ( Figs. [D and jgl), 
we m ade use of the IRAS PSCz catalog (Sau nders et alJ 
l200tf ) , which is the most complete published all-sky galaxy 
redshift catalog. Our analysis essentially consisted of count- 
ing PSCz galaxies within the volume of space in which the 
INTEGRAL or RXTE survey can detect AGN with a given 
hard X-ray (3-20 or 17-60 keV) luminosity, and calculat- 
ing the average galaxy number density in that volume of 
the Universe. To avoid problems with in completeness of 
the P SCz catalog at low Galactic latitudes ijSaunders et alJ 
|2000[ ), we excluded the sky below |6| — 10° from the anal- 
ysis, even though our INTEGRAL AGN sample extends 
down to |6| = 5°. In reality, the situation is somewhat more 
complicated, first because the INTEGRAL and RXTE sur- 
veys are nonuniform, hence the distance out to which a 
source with a given luminosity is detectable depends on the 
direction in the sky, and also because the PSCz catalog it- 
self becomes incomplete at a certain distance that depends 
on galaxy infrared luminosity. Nonetheless, it is not difficult 
to take these effects into account to first approximation. 

Specifically, we first constructed 9 subsamples of IRAS 
PSCz galaxies according to their infrared (60 /im) lumi- 
nosity: LiR > Lcut,,: = 5 X 10^, 10^, 2 X 10^ 5 x 10^, 
10^°, 2 x 10^°, 5 X lO^^ 10", 2 X 10"Lq. These subsam- 
ples are expected to be complete out to Dmax = 23, 32, 
45, 72, 102, 144, 228, 322, 455 Mpc, since the PSCz cata- 
log is highly complete at 60 /xm fluxes higher than 0.6 Jy 
(| Saunders et al .11200 Of ), except near the Galactic plane (we 
ignore here any fc-corrections since we are dealing with low 
redshifts). We then calculated the average number density 
of IRAS PSCz galaxies with Lir > Xcut,i (i = 1.-.9) within 
the volumes of space y(Lhx) probed by INTEGRAL or 
RXTE at different hard X-ray luminosities (17-60 or 3- 
20 keV, respectively): logihx = 41, 41.5, 42, 42.5, 43, 43.5, 
44. For simplicity, in calculating ^(ihx) we assumed power- 
law spectra with a photon index of 1.8. Of all the resulting 
galaxy density estimates we finally took into account only 
those satisfying the condition £>max(icut,i) > DgoiLh^), 
where Dgo is the distance within which 90% of the space 
volume probed by INTEGRAL or RXTE at Lhx is con- 
tained. 



We have thus used a combination of dwarf, normal and 
giant IRAS galaxies to compare the average galaxy num- 
ber densities characteristic of the INTEGRAL and RXTE 
surveys in the broad range of AGN luminosities of our in- 
terest. The result of this comparison is shown in Fig. |4l 
It indicates that the local volume of Universe probed by 
INTEGRAL at L^x ^ 10"^ '^ erg s'^ is a factor of - 1.2 
denser than the corresponding volume probed by RXTE. 
A comparison of the sky exposur e maps of INTEGRAL 
()Krivonos et al.ll2007t ) and RXTE (jRevnivtsev et al.ll2004D 
provides a likely explanation for this: while INTEGRAL 
has spent a lot of time observing such major galaxy concen- 
trations as the Centaurus cluster-Shapley supercluster and 
the Perseus-Pisces supercluster, the RXTE Slew Survey 
has relatively low exposure in these regions. On the other 
hand, we do not find (see Fig. 0]) any significant difference 
in the average galaxy number densities of the INTEGRAL 
and RXTE survey volumes at higher AGN luminosities 
{Lhx ^ lO^"^-^). This is also an expected result, since high- 
luminosity AGN are typically observed by INTEGRAL and 
RXTE at distances ~ 100-200 Mpc and on such spatial 
scales the Universe is believed to be approximately homoge- 
neous. We note however that the limited statistics of distant 
giant galaxies in the IRAS PSCz catalog does not allow us 
to estimate the galaxy number density for the INTEGRAL 
and RXTE surveys at high Lhx to an accuracy better than 
~ 20%. 

To take into account to a first approximation the ap- 
parent effect of the local LSS and assuming that the 
INTEGRAL and RXTE surveys equally over- and un- 
derestimate, respectively, the average density in the local 
Universe, we multiplied the RXTE points of the cumula- 
tive SED of low-luminosity AGN (Fig. (2) by a factor of 
1.1 and divided the INTEGRAL points by the same factor. 
No such corrections were done for the cumulative SED of 
high-hmrinosity AGN (Fig. [3]). 

2.3. Spectral modeling 

We next modeled, using XSPEC ()ArnauJl99l) . the broad- 
band (3-300 keV) cumulative SEDs of local low- and high- 
luminosity AGN obtained in the previous sections by a sum 
of absorbed and unabsorbed power-law spectra with a high- 
energy exponential cutoff: 



E^dN/dE ^ A'^kif{NB..:^E-'^+^ cxp{-E/Ef). 



(4) 



This model is by definition a sum of identical spectra 
absorbed by different column densities Nn i of neutral gas: 

0, 1022-25, 1022-75, 102325, 1023-^5, 1024.25^ 1024-^5 cm-2. 

We fix the weights ki ki = 1) of these columns at the 
relative fractions of AGN with differe nt iVn as measured 
by INTEGRAL (jSazonov et all 120071 ). ignoring the asso- 
ciated uncertainties. The corresponding fractions at low 
(Lhx < 10''3-5 erg s-i) and high (Lhx > 10'*3.5 gj.g g-i) 
minosities are given in TableO We use them accordingly for 
modeling the cumulative SEDs of low- and high-luminosity 
AGN (Fig. [Hand Fig.©. 

The photoabsorption modifiers /(iVu.i) are determined 
using the XSPEC model phabs assuming solar element 
abundances. Note that for substantially Compton thick 
sources (A'h ^ 102"^-^ cm^2) ij^jg simple photoabsorption 
model becomes inadequate, but since the relative fraction 
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Fig. 4. Upper panel: Average number density of IRAS 
PSCz galaxies within the volumes of space probed at |6| > 
10° by INTEGRAL (filled points) and RXTE (open points) 
as a function of AGN hard X-ray luminosity. Each point is 
based on one of 9 subsamples of IRAS PSCz galaxies (from 
top to bottom): Ljr > Lcut = 5 x 10^ 10^ 2 x 10^ 5 x 10^ 
IQi", 2 X IQi", 5 X IQi", 10", 2 x IO^Lq. For a given AGN 
luminosity only one or two best (accurate to within '-^lO- 
30%) density estimates are presented. See text for further 
explanations. Lower panel: Ratio of the densities in the 
INTEGRAL and RXTE volumes as a function Lhx- Also 
shown are the typical ratios, 1.2 and 1.0, at Lhx < 10'*3-5 
(dotted line) and Lhx > 10'*3-5 gj-g g-i (dagjied line), re- 
spectively, which are adopted as correction factors for the 
cumulative SEDs of low- and high-luminosity AGN, respec- 
tively. 



of such sources according to the INTEGRAL survey is small 
(< 10%) we apply this model to all AGN. Also, it would 
be more appropriate to use in equation ^ the true, i.e. 
corrected for selection effects, fractions of AGN with differ- 



Table 3. Observed relative fractions of local AGN 
with different absorption column densities at low and 
high luminosities, u pdated from the values reported by 
ISazonov et al] (|2007D 



logiVH < 22.0 


k (41 < log Lhx < 43.5) 


k (log Lhx > 43.5) 


< 22.0 


0.300 


0.68 


22.0-22.5 


0.175 


0.08 


22.5-23.0 


0.075 


0.04 


23.0-23.5 


0.125 


0.08 


23.5-24.0 


0.150 


0.12 


24.0-24.5 


0.100 


0.00 


> 24.5 
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0.00 
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Fig. 5. Confidence regions for the power-law slope F and 
position of the high-energy cutoff -E/ in the model [equa- 
tion (HI)] for the cumulative SED of low-luminosity (10^^ < 
Lhx < lO'^^-^ erg s-i) AGN. The position of the best-fit 
parameter values is indicated by the cross. 
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Fig. 6. Same as Fig. [51 but for the cumulative SED of high- 
luminosity AGN (Lhx > lO''^-^ erg s^^). 



ent absorption column densities rather than the fractions 
as observed by INTEGRAL in the 17--60 keV energy band 
(Table[3]). However, for the adopted photoabsorption model 
the observed fraction will significantly (by > 10%) under- 
estimate the intrinsic fraction only at A^h ^ 10^^'^ cm~^ 
and since the contribution of such heavily obscured AGN 
to the total SED is expected to be small (<10%) compared 
to the various uncertainties associated with our analysis 
and also because the photoabsorption model itself becomes 
inadequate in the Compton thick regime, we do not try to 
correct our spectral fitting procedure for selection effects. 

The free parameters of our model are the power-law 
slope F, the position of the high-energy cutoff and the 
normalization A, which is allowed to take different values 
for the RXTE and INTEGRAL parts of the spectrum to 
prevent some systematic effects, e.g. associated with the 
cross-calibration of the IBIS/ISGRI and PGA detectors or 
with the local LSS (see i i2.2.ip . from affecting the model- 
ing of the spectral shape. At the same time this allows us 
not to take into account the uncertainties associated with 
the finite sizes of the INTEGRAL and RXTE AGN sam- 
ples [equation ([3])] when fitting the spectrum, since these 
uncertainties are strongly correlated between the different 
energy channels for a given survey and do not affect the 
spectral shape. We do however consider these errors when 
we discuss the final spectral normalization below. 

The above model provides an excellent fit both to the 
cumulative SED of low- luminosity AGN and to that of high- 
luminosity ones. The best-fit parameter values are given in 
Table |4] and the associated confidence regions for the two 
main parameters F and Ef are shown in Figs. [5] and [SI 
The best-fit models are shown in Figs. [D and [31 as well as 
the contributions of unobscured (logA'n < 22), obscured 
(22 < logA'^H < 24) and heavily obscured (logiVn > 24) 
AGN to the total SEDs. 

In summary, the best-fit power-law photon index is 

found to be ^ 1. 7, as expected of AGN (see e.g. 

[Reynolds et"al] 119971 : [Turner et al.lll997[ ). In both studied 
cases, the model prefers to have a high-energy cutoff with 



Ef ~ 100-300 keV, while its detection is only marginal 
(< 2ct) in both cases. Previously, similar values of the cut- 
off energy have been inferred from individual spectra of 
several brightest Sey fert galaxies (in particular NGC 4151, 
[Jourdain et al.[[T993) selected without regard to their lumi- 
nosities or space densities, as well as upon averaging over 
su ch spectra (e.g. Zdziarski et al.[[l995[ : [Perola et al.[[2002[ : 



Mo lina et al][2006ir 



We note the excellent agreement of the derived nor- 
malizations of the low-energy (RXTE) and high-energy 
(INTEGRAL) parts of the cumulative SED of low- 
luminosity AGN, i.e. the whole 3-300 keV spectrum can 
be well fit by a model using a single normalization. Also, 
the ^ 20% difference in the normalizations of the RXTE 
and INTEGRAL parts of the cumulative SED of high- 
luminosity AGN is not statistically significant given the 
uncertainties associated with the photon counting statistics 
and the finite size of the AGN samples (see the first and 
second sets of uncertainties for the A values in Table [IJ . 
The total Icr uncertainties are shown as shaded regions 
in Figs. [2] and [3] and one can see that the INTEGRAL 
and RXTE data are in good agreement with each other. 
Furthermore, part of the apparent mismatch between the 
RXTE and INTEGRAL parts of the cumulative SED of 
high-luminosity AGN may be caused by possible small in- 
homogeneities in the distribution of matter on the 100- 
200 Mpc scales (see ^2XT|l . 

Apart from the relatively simple model given by equa- 
tion ([11), we also fit the measured SEDs by more compli- 
cated models. In particular we tried adding a Compton re- 
flection component. Since this did not lead to a significant 
improvement of the quality of the fit nor did it provide 
interesting constraints on the amplitude of the reflection 
component (due to its strong correlation with the power- 
law slope and the position of the high-energy cutoff), we 
base all of the remaining discussion on our model of a sum 
of absorbed and unabsorbed power-law spectra with a high- 
energy cutoff. 



Table 4. Results of the modeling of the cumulative SEDs of low- and high-luminosity AGN 



erg 


r 


Ef 
keV 


A (INTEGRAL) 

10^** erg 


A (RXTE) 
s-^Mpc"^ 


xVd.o.f. 




1.64 ±0.17 


124+^f 


2.7 ±1.3 ±0.5 


2.7 ±0.6 ±0.6 


3.7/5 




1.74 ±0.08 


258+?^^ 


0.55 ±0.12 ±0.12 


0.44 ± 0.06 ± 0.07 


3.1/5 



2.4. Cumulative SED of the local AGN population 

We have seen that the shapes of the summed SEDs of low- 
and high-luminosity AGN are similar to each other except 
that the spectrum is harder below 20 keV in the former 
case due to the larger fraction of obscured sources at low 
luminosities (Lhx < lO'*'^'^ erg s~^). We can now finally 
determine the SED of the hard X-ray volume emissivity of 
all local AGN with luminosities over 10''^ erg s~^. 

To this end, we first slightly rescale the best-fit spectral 
models presented above. Specifically, we define the ampli- 
tudes of the low- and high-luminosity spectral components 
to be the averages of the corresponding INTEGRAL and 
RXTE values given in Table [H so that 

Alow = 2.7 X 10^* erg s"^ Mpc"^ (5) 
and 

Ahigh = 0.49 X 10^® erg s^^ Mpc"^ (6) 

In Fig. [7] we show by the solid line the sum of our 
low- and high-luminosity best-fit models with the normal- 
izations given by equations ([5]) and ([6]). This is our best 
estimate of the cumulative SED of local AGN. The shaded 
region represents the multifold of spectral models that 
are consistent with the INTEGRAL and RXTE data to 
within Icr and includes an additional uncertainty (±20% 
around the average value) in the spectral normalization, 
which is our estimate of a combination of several uncer- 
tainties, which were already discussed in different sections 
of the paper, namely associated with 1) the finite size of 
the INTEGRAL and RXTE samples of AGN, 2) the local 
LSS (although we have taken it into account to a first ap- 
proximation) and 3) possible incompleteness of both AGN 
samples. The dashed and dotted lines show the contribu- 
tions of the low-luminosity (10^^ < Lhx < lO''^'^ erg s~^) 
and high-luminosity (ihx > lO''"^'^ erg s~^) AGN, respec- 
tively. We note that the latter contribute just ~ 15% to the 
collective SED of local AGN. 

We finally note that blazars are not expected to con- 
tribute more than several per cent to the local hard X- 
ray emissivity at different energies below ~ 200 keV (see 
Fig. [7|). We made this estimate using the same 1/T4iax 
method that we used to derive the cumulative SED of non- 
blazar AGN. Although in making this estimate we did not 
exclude high-redshift flat-spectrum radio quasars present 
in the INTEGRAL and RXTE catalogs, it is essentially 
based on the 2 and 9 nearby (z < 0.1) BL Lacs detected 
by INTEGRAL and RXTE, respectively. We also note that 
the contribution of blazars may well be much more signif- 
icant in gamma-rays {E ^ 200 keV), since emission from 
blazars can extend to very high energies, in contrast to nor- 
mal Seyfert galaxies and quasars. 
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Fig. 7. Cumulative SED of local non-blazar AGN with lu- 
minosities over 10^^ erg s^^. The solid line is our best es- 
timate of this spectrum based on INTEGRAL and RXTE 
data. The shaded region represents the multifold of spectral 
models that are consistent with the data to within la and 
includes the 20% uncertainty in the spectral normalization 
due to a combination of systematic effects discussed in the 
text. The dashed and dotted line show the contributions 
of low- luminosity (10"*^ < Lhx < lO'*'^'^ erg s~^) and high- 
luminosity (Lhx > lO^'^'^ erg s^^) AGN, respectively. The 
two solid curves at the bottom of the plot demonstrate the 
contribution of blazars to the total local emissivity at dif- 
ferent energies, estimated based on INTEGRAL and RXTE 
data, with the distance between the curves corresponding 
to the la uncertainty. 



3. Comparison with the cosmic X-ray background 

In the previous section we derived the cumulative hard X- 
ray SED of local AGN. Specifically, we determined its nor- 
malization, i.e. the volume hard X-ray emissivity of AGN, 
and shape, which was shown to be well described as a sum 
of absorbed and unabsorbed power-law spectra with a high- 
energy cutoff. It is interesting to put these results into the 
broader perspective of the cosmic history of growth of mas- 
sive black holes and the origin of the cosmic X-ray back- 
ground. 



To this end, suppose that the shape of the cumulative 
AGN SED does not evolve with redshift, which implies a 
constant relative contribution of obscured sources, while 
its normalization, i.e. the luminosity density of AGN, does 
experience evolution. The collective emission of AGN from 
all redshifts observed at z = will then have the s pectrum 
(for a flat cosmoslogy, see e.g. ISazonov et al.|[200^ 



e{z)S{{l + z)E) 



dz. 



(7) 



where S{E) is the cumulative SED of local AGN (Fig. E]) 
and the function e(z) describes the evolution of the AGN 

luminosity de nsity. 

Following ISazonov et al.l (|2007D . we assume that the 
hard X-ray luminosity density of AGN evolves similarly to 
the luminosity density of AGN with Lhx ^ 10'*^ erg s~^ 
in the rest-f r ame 2-8 keV energy band as reported by 
iBarger et al.l (|2005D . specifically that e(z) is bound between 
two limiting functions: 



ei(z) cx 



and 



€2(2) cx 



{l + zf-^ 



(H-z)^■^ 
^2(1), 



z < 1 
z > 1 



z < 1 
z > 1. 



(8) 



(9) 



This formulation reflects the significant uncertainty in our 
knowledge of the AGN evolution at z ^ 1, whereas the 
rapid, approximately power-law evolutio n between z = 
and z ~ 1 is established relatively well (lUeda et al.l l2003t 
IBarger et al]|2005l: iLa Franca et al.ll2005D . 

We now use equation ([7]) to convolve our locally deter- 
mined composite AGN SED (Fig. [7]) with the two limiting 
evolution functions ei and €2 (up to Zmax = 6). The re- 
sulting SEDs, representing the collective emission of AGN 
from all redshifts, are shown in Fig. [Sjand Fig. [51 respec- 
tively. The shaded uncertainty regions in these figures have 
been derived by propagating the uncertainties associated 
with the local AGN SED shown in Fig. [71 including the la 
uncertainty in the spectral shape and the estimated 20% 
uncertainty in the volume cmissivity of local AGN. 

We see that within the fairly large uncertainties the 
shape of the predicted redshift-integrated AGN SED is in 
good agreement with that of the CXB for both limiting 
scenarios of AGN evolution at z > 1 and thus for the true 
evolution constrained between ei(z) and e2(z). As concerns 
the normalization, depending on the character of evolution 
at z > 1 our model is between being consistent with the 
observed CXB intensity (in the case of e2{z)) and under- 
estimating it by some 30% (in the case of ei(z)). We also 
note that by adopting the functional form e(z) = (1 -I- z)^^^ 
at z < 1 we essentially fixed the ratio of the AGN lumi- 
nosities densities at z = 1 and z = at 9.2. In reality , 
this ratio is known with accuracy ~ 20%lUeda et~aIll2003L 
IBarger et all 120051: iLa Franca et all |2005| and this uncer- 
tainty of course translates into a similar uncertainty in 
our predicted spectrum. The calculated SEDs are com- 
pared in Fig. |5| and ID with the CXB spectrum measured 
by INTEGRAL (IChurazov et"al]l2007D . 

We note that if the simple scenario of evolution con- 
sidered here is correct, we may expect the calculated SED 
to somewhat understimate the CXB intensity. Indeed, the 
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Fig. 8. Result of convolution of the cumulative SED of local 
AGN (Fig. [7]) with the redshift evolution function of the 
AGN luminosity density described by equation ([5]). The 
solid line is our best estimate of this redshift-integrated 
AGN SED. The shaded region represents the Icr uncertainty 
in the spectral shape combined with the 20% uncertainty 
in the normalization. The points with error bars show the 
CXB spectrum measured with the JET-X, IBIS/ISGRI a nd 
SPI instruments on INTEGRAL (jChurazov et al.ll2007f ). 



calculation above rested on our estimate of the local AGN 
emissivity, which was based on non-blazar AGN with lu- 
minosities higher than 10^^ erg s~^. An additional signifi- 
cant contribution to the local emissivity is expected to come 
from AGN with st ill lower luminosities (Lhx < 10'*"'^ ergs""'^) 
(lElvis et al.l I1984D and from blazars (see ^2A\ above). In 
fact, we recently used RXTE observations to study the cor- 
relation of large-scale variations of the CXB intensity with 
the spatial distribution of galaxies in the local Universe, 
which allowed us to estimate the combined emissivity of 
all X-ray sources (AGN, star-forming and normal galax- 
ies, and clusters of galaxies etc) in the local Universe at 
9 ±4) X 10^^ erg s~^ M pc~^ in the 2-10 keV energy band 
Revnivtsev et al.ll2008t) . On the other hand, by integrating 
the luminosity function of nearby AGN measured by RXTE 
(Sazonov & Revnivtsev 2004) and applying the corrections 
described in H2.2[ we find that the combined emissivity of 
AGN with Lhx > lO''^ erg s'^ is ~ 6 x lO^^ erg s"! Mpc'^ 
Therefore, low-luminosity (Lhx < 10**^ erg s ^^) X-ray 
sources may indeed contribute ~ 30% or even more to the 
total local X-ray emissivity (at least in the 2-10 keV energy 
band and probably somewhat less at higher energies taking 
into account that clusters of galaxies and normal galaxies 
produce softer X-ray emission than AGN). 

We also point out that our calculation was based 
on the assumption that the evolution of the total lu- 
minosity density of AGN follows the observed evolution 
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Fig. 9. Same as Fig. [9l but for the scenario of evolution 
given by equation 



of the 2-8 ke V luminosity density of AGN with ihx 
10*^ erg s~^ (|Barger et alj [2005 1 . The inclusion of lower- 
luminosity (I/hx ^ 10^^ erg s~^) AGN, which are not yet 
accessible to observation at z ^ 1 and whose abundance 
is thus poorly known, could also modify the result of our 
calculation. 

We finally note the well-known fact that if the luminos- 
ity density of AGN were the same at high redshifts as it is 
at z = 0, only ^ 20% of the CXB could be explained. 



4. Discussion and conclusions 

In this paper we used the INTEGRAL and RXTE all-sky 
hard X-ray surveys to calculate the spectral energy distri- 
bution of the collective emission of local AGN in the broad 
range 3-300 keV, properly taking into account the relative 
contributions of AGN with different luminosities and ab- 
sorption column densities. 

We first performed stacking spectral analyses separately 
for AGN with low (lO^^ < Lhx < lO''^-^ erg s'^) and high 
(Lhx > 10^^'^ erg s^^) luminosities and found that both 
stacked SEDs are consistent with having the same shape at 
energies above 20 keV - in particular, a high-energy cutoff 
is marginally detected above ~100-200 keV, whereas the 
cumulative SED of low-luminosity AGN is harder at ener- 
gies below 20 keV due to the larger fraction of obscured 
AGN compared to higher luminosities. We then summed 
up our best-fit models for the stacked spectra of low- and 
high-luminosity AGN to obtain the collective SED of all lo- 
cal non-blazar AGN with Lhx > 10'*^ erg and estimated 
the associated statistical and systematic uncertainties. 

As a first attempt to apply the derived cumulative SED 
of local AGN to studying the history of black hole growth 
in the Universe, we demonstrated that this SED is consis- 



tent with the cosmic X-ray background spectrum if both 
the spectral shape of the collective hard X-ray emission 
of AGN and the relative fraction of obscured AGN remain 
constant with redshift while the total AGN luminosity den- 
sity undergoes strong evolution between z ^ 1 and z — 0. 

This simple model underpredicts the observed CXB in- 
tensity by ~ 0-30% depending on the loosely constrained 
evolution of AGN at z > 1 (and to a lesser degree on the 
evolution at z < 1). As we discussed in 321 the missing 
CXB flux can also possibly be attributed to blazars and 
to low-luminosity AGN (Lhx ^ 10^^ erg s~^ at z = and 
Lhx ^ 10^^ erg s^^ at z ^ 1~2), not included in our model, 
since such objects are known to be abundant and ener- 
getically important at least in the local Universe. In this 
connection we note that there may be an increased fraction 
of Compton thick (iVn ^ 10^'* cm~^) objects among such 
low-luminosity AGN compared to the ~ 15% fraction at 
ihx ^ 10^^ erg s-i observed by INTEGRAL. Indeed, some 
studies suggest that the intrinsic fraction of Compton thick 
AGN in the local Universe may be as hig h as ^ 50% (e.g. 
iRisaliti e t al. 19991: iGuainazzi et al.ll2005[ ) and in the most 
obscured of these objects even the hard X-ray flux will be 
strongly depressed. 

Turning the above argument around, we can conclude 
that AGN with Lhx ^ iC^ erg s^^, which are undetectable 
even in the deepest extragalactic surveys with Chandra and 
XMM-Newton, cannot provide a dominant contribution to 
the X-ray emissivity of the Universe at redshifts z ~ 1-2, 
otherwise they would produce too much X-ray background. 

How does this assumed scenario of AGN evolution com- 
pare with observational data obtained by X-ray surveys at 
energies below ~ 8 keV? First of all, a number of studies 
have shown that AGN undergo approxi mately pure lumi- 
nos ity evolution be tween z ~ 1 and z = (jUeda et al.ll2003t 
B arger et al.ll2005[ ). which means that the AGN luminosity 
function has shifted to lower luminosities by nearly an or- 
der of magnitude. More fundamentally, it has been demon- 
strated that black hole growth and galaxy formation move 
steadily and in parall el to low er mass scales since z 2 up 
to z = (e.g. Heckman et al.|[2004i V 

This picture of "cosmic downsizing" or " antihierarchi- 
cal evolution" is consistent with the scenario considered in 
this paper and in fact in our calculations we adopted the 
functional form of the AGN luminosity density evolution 
between z = and z ~ 1 from the recent X-ray surveys. As 
concerns the relative contribution of obscured AGN to the 
cumulative hard X-ray emissivity at a given redshift, which 
is constant in our model, observations seem to indicate that 
there is indeed no signific ant evolution of t his fraction be- 
tween z = and z ~ 1 (jUeda et al.l 120031 ). Therefore, all 
the available AGN observations at z ^ 1 appear to well fit 
in the simple scenario of evolution considered in this paper. 

However, X-ray observations also suggest that the char- 
acter of AGN evolution changes at higher redshifts (z ^ 1). 
Namely, the observed AGN evolution is much better de- 
scribed by a luminosity dependent densit y evolution model 
than by a pure l uminosity evolution one ^Ueda et al.ll2003l : 
iLa Franca et all [2005: GiUi et al.|[2007t ). Also the relative 
contribution of obscured sources to the total AGN lumi- 
nosity density may have undergone substantial evolution 
at z ^ 1, although this remains a very controversial is- 



r y controv 

sue (La Fra nca et al.ll2005tlAkvlas et al.ll2006l : iTozzi et all 
2006; Treister fc Urrvll2006( ). Although the different char- 
acter of AGN evolution at high redshifts with respect to 



z < 1 suggests that our simple model needs to be modified, 
such a revision would remain almost unnoticed in view of 
the fairly large uncertainties in the cumulative SED of local 
AGN. Indeed, as has been shown in Sj3l the current uncer- 
tainty in the AGN evolution at z > 1 leads to less than 
~ 30% uncertainty in the predicted CXB spectrum, which 
is of the same order as the uncertainties in the cumulative 
SED of local AGN obtained from INTEGRAL and RXTE 
data. 

We also note that in reality the hard X-ray spectra 
of distant quasars may somewhat differ from those of the 
nearby Seyfert galaxies that make up our local cumulative 
SED, since, although AGN spectra are not expected to di- 
rectly depend on cosmological redshift, they can depend 
on physical parameters of the accretion disk around the 
central massive black hole, determined by the black hole 
mass, accretion rate (Shakura & Sunyaev 1976) and spin. 
Observations with future hard X-ray telescopes will permit 
direct tests of whether the high-energy spectra of quasars 
are similar to those of local AGN or not. 

The main result of this work is that for the first time 
a direct comparison has been made between the cumula- 
tive hard X-ray SED of the local AGN population and the 
CXB spectrum, which demonstrated that the commonly 
accepted paradigm of the CXB being a superposition of 
AGN is generally correct. Improved measurements of the 
cumulative SED and evolution of AGN by current and fu- 
ture X-ray and hard X-ray astronomy missions will make it 
possible to obtain tighter constraints on the cosmic history 
of black hole growth and the AGN unification paradigm. 
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